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Philosophy  of  Theoretical  Research  In  Turbulence 


The  objective  of  this  research  ef  ort  Is  to  construct  a  purely  theoretical 
foundation  for  the  phenomenon  of  turbulent  fluid  flow  and  heat  transfer.  The 
starting  point  In  this  effort  Is  the  Principal  Investigator's  view  that,  as 
an  unexplained  physical  phenomenon,  turbulence  represents  a  new  physics,  and 
that  a  new  phsycis  demands  a  new  theory  which  must  be  judged  based  on  entirely 
new  experiments.  This  view  recognizes  that,  historically,  the  classical 
Navter-Stokes  equations  have  been  shown  to  be  Inadequate  as  a  theory  of  tur¬ 
bulence.  One  hundred  years  of  consistent  inability  to  avoid  empiricism  while 
proposing  "theories"  with  the  Navier-Stokes  equations  as  starting  point,  is 
the  best  indicator  that  an  entirely  new  theoretical  viewpoint  Is  needed. 

Indeed,  based  on  Its  record  of  accomplishments,  the  theory  embodied  in  the 
classical  Navier-Stokes  equations*  emerges  as  only  a  theory  of  laminar  flow. 

The  fundamental  theoretical  question 

For  a  purely  theoretical  viewpoint  to  be  successful  In  explaining  and 
predicting  the  occurence  of  "turbulence".  It  must  be  able  to  account  for  the 
most  basic  ,-snd  elementary  feature  of  turbulent  flow.  This  feature  is  un¬ 
questionably  the  "eddy",  l.e.,  that  finite-size  portion  of  the  flow  field 
which  rotates  and  -  as  a  unit  -  Hows  the  bulk  motion  of  the  field.  The 
fundamental,  zeroth-order  question  to  be  answered  Is  then: 


What  "makes"  the  FIRST  EDDY,  in  a  field  which 
originally  contains  absolutely  no  sign  (trace) 
of  turbulent  motion? 


*The  principle  of  mechanical  equilibrium  of  a  point-size  fluid  packet  in 
communication  with  neighboring  point-size  packets. 
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It  Is  worth  pointing  out  that  this  fundamental  question  Is  totally  dis¬ 
regarded  by  "theories"  whtch  proceed  from  the  time-averaged  Navter-Stokes 
equations,  a  procedure  pioneered  by  Osborne  Reynolds  himself  100  years  ago. 

Such  analyses  take  the  presence  (birth)  of  eddies  for  granted,  hence,  they 
are  by  definition  empirical  (empiricism  Is  the  thinking  process  which  has  as 
starting  point  an  unexplained  physical  observation).  It  is  for  this  reason 
that  the  current  search  for  "generally  applicable  models"  for  time-averaged 
eddy  transport  Is  destined  to  fail:  to  search  for  such  models  Is  the  philo¬ 
sophical  equivalent  of  claiming  to  understand  the  behavior  of  a  society  while 
negating  (disregarding  the  importance  of  understanding)  the  Individual. 

The  fundamental  question  formulated  above  appears  to  have  been  avoided 
by  turbulence  theorists  over  the  years.  Of  course,  there  are  explanations 
for  this  unfortunate  record.  For  example.  It  Is  generally  acknowledged  that 
the  interface  between  the  understood  (laminar)  and  the  not  understood  (turbulent) 
behavior  of  fluid  flow  is  described  satisfactorily  by  the  Theory  of  Hydro- 
dynamic  Stability.  Yet,  the  starting  point  in  any  hydrodynamic  stability 
analysis  is  the  postulate  that  "disturbances"  of  known  size  exist,  super¬ 
imposed  on  the  "undisturbed"  (laminar)  field.  With  regard  to  explaining 
the  origin  of  turbulence,  the  "theory"  of  Hydrodynamic  Stahl lity  is  nothing 
but  an  empirical  exercise*,  the  game  of  watching  what  happens  (In  time)  to 
turbulence  which  is  already  present.  Hydrodynamic  stability  studies  negate 
the  importance  of  answering  theoretically  the  fundamental  question  concerning 
the  natural  origin  of  disturbances..  In  '.his  way,  the  field  of  Hydrodynamic 
Stability  illustrates  the  Intrinsic  inability  of  the  Navler-Stokes  theory  to 
explain  turbulence. 

*An  experiment  which,  historically,  was  first  pursued  analytically  and,  more 
recently,  numerically  and  in  the  laboratory  (e.g.,  the  loudspeaker  harmonic 
forcing  of  an  "undisturbed"  flow). 
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There  are  other,  equally  Important  reasons  for  why  the  fundamental 
question  of  theoretical  research  in  turbulence  has  not  been  asked.  Ore 
reason  Is  the  fact  that,  traditionally,  old  theories  and  old  theoreticians 
hide  their  Inability  in  the  face  of  a  new  physics  by  inventing  and  dissemina¬ 
ting  a  "new"  terminology  which  is  misleading  and  self-serving.  Note  the 
postulate  of  "disturbances  In  a  laminar  field"  employed  In  Hydrodynamic 
Stability  studies  which  go  on  to  claim  that  they  shed  light  on  the  "transition" 
from  laminar  to  turbulent  flow.  In  reality,  the  "transition"  has  already 
taken  place,  before  the  start  of  the  stability  experiment.  The  initial 
condition,  interestingly  described  as  "disturbances  superimposed  on  a  laminar 
flow  field",  is  nothing  but  "unexplained  turbulence".  The  lack  of  emphasis 
science  educators  place  on  seeing  through  the  language  coined  by  any  theory, 
oid  and  new,  is  the  major  contributor  to  lack  of  progress  in  theoretical  re¬ 
search,  hence,  to  the  conversion  of  science  to  religion. 

The  elasticity  of  turbulent  streams. 

Turbulence  Is  a  characteristic  behavior  of  Inviscld  flow.  In  order  to 
see  this,  consider  the  simple  experiment  shown  in  the  figure.  A  stream  of 
fixed  flowrate  m  flows  through  a  pipe  of  fixed  diameter  D.  If,  chemically, 
the  experimentalist  always  decreases  the  fluid  viscosity  y,  the  flow  will 
eventually  shift  from  the  predictable  laminar  flow  to  the  unpredictable 
turbulent  behavior.  At  the  same  time,  it  can  be  said  that  the  flow  shifts 
from  viscid  to  the  inviscld  flow  limit.  To  be  sure,  nothing  is  viscid  (or 
inviscld)  everywhere  and  forever;  however,  it  Is  clear  that  the  viscidity  of 
the  stream  as  a  whole  (as  a  finite  size  system)  decreases  as  y  decreases. 


experiment  • 

VISCID 

1 

y  decreases  j  INVISCID 

laminar 

I  turbulent 

'  high  Re 

|  zero  entropy  generation 

|  elastic  (conservative) 

Another  important  aspect  of  the  invlscld  flow  limit,  the  limit  of  tur¬ 
bulent  behavior,  is  the  fact  that  the  entropy  generation  rate  steadily  de¬ 
creases  as  y  decreases.  The  pipe  becomes  a  thermodynamic  system  filled  with 
conservative  material  incapable  of  irreversibility.  Such  a  system  possesses 
ilastic  properties.  Indeed,  if  the  pipe  used  in  the  experiment  Is  highly 
flexible,  then  it  is  possible  to  store  work  in  the  system  (the  pipe)  by 
temporarily  curving  it.  Just  like  a  garden  hose  delivering  a  high  (turbulent) 
flowrate,  the  elastic  stream-pipe  system  has  the  natural  tendency  to  spring 
back. 

The  thinking  experiment  illustrated  in  the  figure  leads  to  a  fundamental 
conclusion  in  the  search  for  the  origin  of  turbulence: 

Turbulence  is  the  behavior  of 
invlscld  (elastic)  streams. 


-  H 
e.  .J 
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The  sinusoidal  shape  of  invlscid  (elastic)  streams. 

Wa  learn  more  about  turbulence,  as  the  natural  behavior  of  Invite  Id 
streams*,  by  taking  an  unbiased  look  at  the  tendencies  of  highly  tnvtsclu 
streams  we  see  in  nature.  The  first  piece  to  look  Is  the  flow  of  rivers 
where  we  recognize  umeandersM  as  a  characteristic  of  rivers  of  all  sizes. 
Meanders  represent  a  highly  regular  geometric  pattern  resembling  a  sinusoid 
with  exaggerated  elbows.  But,  most  Importantly,  there  appears  to  be  a  uni¬ 
versal  proportionality  between  the  wavelength  of  the  sinusoidal  course  and 
the  width  of  the  river. 

Meanders  are  a  characteristic  of  all  streams  in  the  turbulent  flow  limit. 
Just  think  of  the  snake-1  Ike  shape  of  all  turbulent  plumes,  the  "buckled" 
shape  of  lightning  and  the  meanders  of  rivulets  (water  trick) Ings  on  a  shower 
curtain,  car  windshield,  etc.).  We  conclude  that: 


Natural  inviscid  streams  which 
are  not  constrained  by  rigid 
walls  have  the  property  to 
meander,  so  that 
the  meander  wavelength 

th.  width  '  Ul’,v,r"1  constant 


*aTso  cal  led  "high  Reynolds  number"  flows. 
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The  buckling  property  of  inviscld  streams. 

To  summarize  the  preceding  discussion,  we  found  that  finite-size  in- 
viscid  streams  are  elastic  (springy)  and  that,  if  unconstrained,  inviscld 
streams  assume  a  regular,  sinusoidal  shape.  If  we  think  only  about  the 
finite-size  space  (column)  occupied  by  the  inviscld  stream,  and  not  about 
the  fluid,  the  above  conclusions  can  be  restated  as  follows: 

The  column  contains  springy  material  and. 

If  free  to  deform  laterally,  assumes  a 
regular  sinusoidal  shape. 

This  is  precisely  the  physics, i  observation  which  leads  to  the  Buckling 
Theory  of  slender  columns  in  the  field  of  Strength  of  Materials.  Indeed, 
the  stream  column  possesses  the  previously  unknown  property  of  buckling,  as 
shown  in  the  remaining  sections  of  this  report. 

The  newly  discovered  buckling  property  of  inviscld  streams  (fluid  fibers 
of  finite  thickness)  accounts  for  the  "turbulent"  motion  of  fluids.  First 
examples  of  how  this  property  explains  turbulence  are  given  in  the  1st  and 
3rd  papers  in  this  report.  From  the  outset,  It  is  important  to  recognize 
that  the  buckling  property  of  inviscld  flow  Is  not  another  "model"  of  the  kind 
proposed  frequently  in  contemporary  turbulence  research.  The  buckling 
property  of  inviscld  streams  is  a  purely  theoretical  discovery,  an  integral 
part  of  newtonian  fluid  mechanics  which  has  eluded  researchers  until  now.  A 
reexamination  of  the  1st  paper  which  develops  the  buckling  theory  of  inviscld 


fibers  of  finite  thickness,  reveals  that  the  buckling  property  is  an  im¬ 
mediate  consequence  of  the  classical  Bernoulli  equations.  The  new  feature 
of  the  Buckling  Theory  is  the  fact  that  the  system  of  interest  is  of  finite 
size.  This  represents  a  drastic  departure  from  previous  theoretical  research 
in  fluid  mechanics  which,  starting  with  the  Nav ter-Stokes  equations,  focuses 
on  the  equilibrium  of  a  point-size  system. 

Built  into  the  Buckling  Theory  is  the  new  philosophy  that  a  finite-size 
portion  of  the  flow  possesses  a  global  property  (buckling).  Consequently, 
the  fluid  contained  in  a  fiber  evolves  as  a  finite-size  entity,  not  as  a 
point  size  fluid  packet  in  laminar  flow.  Inviscid  flow  is  a  fibrous  flow: 
fluid  fibers  are  distinct  portions  of  the  flow  field,  defined  by  surfaces  of 
relative  motion  (e.a.,  jets,  wakes,  shear  layers).  Fluid  fibers  buckle  and 
the  elbows  of  the  buckled  shape  mix  with  the  ambient  giving  birth  to  eddies. 

The  new  eddies  thicken  the  old  fiber  Into  a  new  fiber  which  also  buckles. 
Published  Papers. 

No.  1.  A.  Bejan,  "On  the  Buckling  Property  of  Inviscid  Jets  and  the  Origin 

of  Turbulence."  Letters  in  Heat  and  Mass  Transfer,  Vol .  8,  pp.  187-194, 
May-June  1981. 

No.  2.  A.  Bejan,  "Comments  on  Viscous  Buckling  of  Thin  Fluid  Layers," 

Physics  of  Fluids,  Vol.  24,  No.  9,  pp.  1764,1765,  September  1981. 

No.  3-  A.  Bejan,  "Theory  of  Instantaneous  Sinuous  Structure  in  Turbulent 
Buoyant  Plumes,"  Warme-und  Stoffubertragung,  Vol.  1 6  ,  March  1982 
(to  appear) . 
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Paper  no.  1 


ON  THE  BUCKLING  PROPERTY  OF  INVISCID  JETS  AND 
THE  ORIGIN  OF  TU13ULENCE* 

Adrian  Bejan 

Department  of  Mechanical  Engineering 
University  of  Colorado,  Boulder,  Colorado 

Abstract 

This  paper  outlines  the  analogy  which  exists  between  inviscid  jets 
and  elastic  columns  in  axial  compression.  It  is  shown  that  straight  in¬ 
viscid  jet  columns  possess  the  property  of  sinusoidal  infinitesimal 
buckling.  The  buckling  wavelength  scales  with  the  transversal  dimension 
of  ths  jvjt.  The  repeated  buckling  and  breakup  of  the  jet  column  is 
responsible  for  the  observed  whiplash  motion  of  turbulent  jets.  The 
buckling  theory  predicts  correctly  the  natural  frequency  of  the  whiplash 
motion  and  the  Reynolds  number  for  the  laminar-turbulent  transition  in 
free  jet  flow. 


* 

published  in  Letters  in  Heat  and  Mass  Transfer,  Vol.  8,  pp.  187-194,  1981. 
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One  aspect  of  jet  turbulence  receiving  increased  attention  is  the 

large  scale  "orderly"  structure  with  a  length  scale  of  the  same  order  as 

the  jet  diameter.  Crow  and  Champagne  [l]  showed  that  as  the  Reynolds 

2  3 

number  increases  from  10  to  10  ,  the  whiplash  motion  of  jets  evolves 
from  a  sinusoid  to  a  helix  and  eventually  to  a  train  of  axisymmetric 
waves,  similar  observations  have  been  reported  by  Reynolds  [2],  The 
characteristic  snake-like  shape  of  a  turbulent  round  jet  is  shown  very 
clearly  in  Fig.  1. 

The  purpose  of  this  paper  is  to  offer  a  theoretical  explanation  for 
the  observed  large-scale  periodic  structure  of  turbulent  jets.  The  ex¬ 
planation  is  founded  on  a  very  interesting  analogy  which  exists  between 
jet  flows  and  slender  elastic  columns  in  axial  compression.  It  may  be 
recalled  that  Euler's  theory  of  infinitesimal  buckling  (indifferent 
equilibrium)  in  axially  compressed  columns  rests  on  only  two  premises  [4-6] ; 
(i)  the  slender  column  is  straight  and  in  axial  compression; 

(ii)  if  subjected  to  a  separate  bending  test,  the  column  develops 
in  its  cross-section  a  resistive  bending  moment  which  is  pro¬ 
portional  to  the  induced  curvature. 

slender  elastic  columns,  of  course,  meet  these  two  conditions.  However,  it 
is  shown  in  the  next  section  that  exactly  the  same  conditions  are  met  by 
inviscid  streams  discharging  freely  into  larger  reservoirs.  Consequently, 
the  column  (control  volume)  occupied  by  an  inviscid  stream  buckles  sinu¬ 
soidally,  and  the  stream  mixes  periodically  with  the  stagnant  ambient. 

From  this  result,  we  conclude  that  the  natural  property  of  inviscid  streams 
is  to  follow  a  sinusoidal  (meandering)  path  as  they  travel  through  a 
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stagnant  ambient.  This  natural  property  of  inviscld  flow  is  the  basis  for 
the  unexplained  turbulent  behavior  of  fluids. 

The  static  Equilibrium  of  the  Jet  Envelope 

Consider  a  straight  inviscid  jet  of  density  p  ,  uniform  velocity 
V  and  cross-sectional  area  A  ,  as  shown  in  Fig.  2.  The  static  pressure 
inside  the  jet  and  in  the  ambient  fluid  is  Pq  .  Imagine  now  a  stationary 
envelope  (control  surface)  which  surrounds  a  certain  length  of  the  jet. 

This  envelope  and  two  transveral  end-cuts  define  a  stationary  control 
volume . 

In  the  spirit  of  the  thermodynamics  of  open  flow  systems  [7],  the  only 
forces  which  act  on  this  control  volume  are  the  inlet  and  outlet  compres¬ 
sive  forces 

C  =  p  A  V2  (1) 

Forces  C  are  shown  schematically  in  the  lower  half  of  Fig.  2  where  the 
control  volume  is  symbolized  by  the  solid  line.  At  this  point  we  conclude 
that  the  fluid-filled  column  represented  by  the  control  volume  satisfies 
condition  (i)  for  buckling. 

when  the  control  volume  is  slightly  curved,  each  face  of  the  trans¬ 
versal  cut  is  exposed  not  only  to  a  compressive  force  but  also  to  a  bending 
moment.  Consider  the  separate  bending  test  in  which  the  jet  is  held 
(forced  to  flow)  in  a  slightly  curved  duct.  The  radius  of  the  curvature 
of  the  duct,  Rw  ,  is  infinitely  greater  than  the  transversal  dimension  of 
the  jet,  D.  Bernoulli's  equation  for  a  streamline  [8]  dictates 

r  P  V2  +  P  -  -  p  v2  (z)  +  P(z) 

2  O  2 


(2) 
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where  the  right-hand-side  of  the  equation  corresponds  to  the  curved  duct. 
Coordinate  z  is  measured  radially  from  the  jat  centerline  towards  the 
center  of  curvature.  Radial  equilibrium  of  the  jet  fluid  in  the  curved 
duct  requires  also  [9] 

ov2(z)  3P 

-nr~“i7  (3) 

oo 

In  the  limit  of  vanishingly  small  curvature,  D/R  -+  0  ,  Equations  (2)  and 

00 

(3)  yield 


v  (z)  =  V  ( 1  +  ^-> 

OO 


(4) 


P(z) 


p  -  Pv2z 

°  R» 


(5) 


The  bending  moment  acting  over  the  cross-section  is 


M  = 


(pv2  +  P)  z  dA 


pV2I 

R 

00 


(6) 


where  I  = 


z2dA  is  the  area  moment  of  inertia  of  the  jet  cross-section. 


In  conclusion,  the  cross-sectional  bending  moment  M  of  a  nearly 

straight  jet  column  is  proportional  to  the  curvature  of  the  column,  1/R  . 

00 

This  means  that  l*:e  inviscid  jet  column  satisfies  condition  (ii)  for  the 
infinitesimal  buckling  of  a  straight  column. 

Consider  now  the  static  equilibrium  of  the  straight  jet  envelope 
shown  in  Fig.  2.  Clearly,  the  axial  compressive  forces  C  balance  each 
other.  In  addition,  the  excentricity  bending  moment  CY  —  however 
small  —  must  be  balanced  at  all  times  by  the  cross-sectional  bending 
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moment  [10] 


CY  +  M  -  M (x)  =  0 
o 


(7) 


Noting  that  1/R^  =  -d2Y/dx2  ,  the  rotational  equilibrium  condition  (7) 
constitutes  a  differential  equation  for  the  equilibrium  centerline  Y(x)  , 


pV2I 


d2Y 

dx2 


+  pV2AY  +  M 

o 


0 


(8) 


0  at  x  -  0  is 

(9) 

where  the  amplitude  K  is  indeterminate  and  infinitely  small  compared  with 
the  transversal  length  scale  /i/A  or  D  . 

We  arrive  at  the  important  conclusion  that  the  equilibrium  centerline 
of  the  jet  envelope  is  a  sinusoid  of  infinitely  small  amplitude.  The 
wavelength  of  this  trajectory  scales  only  with  the  transversal  dimension 
of  the  jet, 

A  *  2tt/T7a  (10) 


The  solution  satisfying  the  nozzle  conditions  Y  =  Y'  = 
Y (x)  =  k[cos (x/a/I)  -  lj 


For  example,  A/D  =  n/2  for  a  round  jet  of  diameter  D  ,  and  A/D  =  n//3 
for  a  flat  (two-dimensional)  jet  of  thickness  D  .  An  analysis  of  the 
meandering  contour  of  the  round  jet  of  Fig.  1  yields  A/D  =  1.2  which  is 
in  good  agreement  with  the  theoretical  value  of  1.57  (it  is  very  likely 
that  the  24%  discrepancy  between  the  two  values  is  caused  by  the  fact  that 
the  counter  shown  in  Fig.  1  overestimates  the  real  diameter  of  the  turbu¬ 


lent  jet) . 
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Formation  of  Large-Scale  Vortices 


The  analogy  between  the  infinitesimal  buckling  of  elastic  columns 
and  inviscid  jet  envelopes  terminates  with  the  equilibrium  centerline 
given  by  Equation  (9) .  Whereas  in  slender  elastic  columns  the  equilibrium 
of  the  small-amplitude  sinusoid  is  indifferent  [10] ,  in  an  inviscid 
jet  the  equilibrium  is  unstable.  The  unstable  equilibrium  of  the  inviscid 
column  is  described  very  well  by  the  classical  theory  of  hydrodynamic 
instability  [11].  The  slightest  deviation  of  the  jet  from  its  rectilinear 
shape  leads  to  the  formation  of  lateral  "lift"  forces  which  consistently 
tend  to  amplify  the  deformation.  As  a  result,  the  inviscid  jet  breaks 
up  periodically  as  its  overextended  elbows  penetrate  and  mix  with  the 
stagnant  ambient.  The  degenerated  elbow  region  becomes  a  "large  scale 
turbulent  structure"  [1]  which  continues  to  move  downstream  with  a  speed 
of  order  V/2  . 

The  repeated  buckling  and  breakup  phenomenon  accounts  for  the  observed 

whiplash  and  fluctuating  motion  of  turbulent  jets.  The  period  of  this 

fluctuation  scales  with  the  buckling  time  t  *  A/V  ,  i.e.,  with  the  time 

B 

of  fluid  travel  between  successive  elbows  (breaks)  in  the  jet  column. 

The  natural  frequency  of  jet  fluctuation  can  be  expressed  in  dimensionless 
form  as  a  Strouhal  number 


D  _  D 
t  V  “  A 


(ID 


which,  as  shown  in  the  preceeding  section,  has  a  *•  lue  of  order  0.5. 

This  order  of  magnitude  estimate  agrees  very  well  with  experimental  measure¬ 
ments  of  the  natural  frequency  cr  turbulent  jets  exposed  to  a  range  of 
external  excitation  frequencies.  For  example,  Bechert  and  Pfizenmaier  [12] 
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reported  St  *  0.5  for  maximum  amplification  of  broadband  jet  noise. 
Most  recently#  Acton  [ J 3 J  conducted  a  computational  simulation  of  a 
round  turbulent  jet  by  modeling  the  shear  layer  as  a  succession  of  dis¬ 
crete  vortex  ring  elements.  Acton  showed  that  the  natural  periodicity 
of  the  round  jet  was  St  *  0.47  ,  and  that  the  jet  was  most  sensitive 
when  forced  at  approximately  the  same  frequency  (st  «  0.5). 


The  Transition  to  Turbulence 


The  buckling  property  and  fluctuating  nature  of  inviscid  jets  pro¬ 
vide  a  theoretical  ba?is  for  predicting  the  transition  to  turbulence. 
The  transition  from  laminar  jet  flow  to  turbulent  flow  occurs  when  the 
statxcnary  ambient  is  no  longer  capable  of  viscously  communicating  with 
the  jet.  The  viscous  communication  time  between  the  jet-ambient  inter¬ 
face  and  the  jeu  centerline  follows  from  Stokes’  first  problem  1.14) 


t  » 
v 


d_L 

16v 


(12) 


The  jet  is  free  to  buckle,  i.e.,  to  get  out  of  hand,  rf  the  viscous  diffu¬ 
sion  time  tv  is  longer  than  the  buckling  and  breakup  time  t  .  Defining 
the  buckling  number  N  as  the  ratio  of  these  two  characteristic  times 

O 

of  the  f . ow  configuration, 


N 


B 


(13) 


we  have  a  criterion  which  predicts  the  fluctuating  (turbulent)  behavior 

of  the  jet.  Using  Equation  (12)  and  t  *  1/V  ,  the  buckling  number  can 

B 


he  written  also  as 
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N  .  VD  D^A 
B  v  16 


(14) 


For  a  round  jet,  the  transition  criterion  N  >  l  becomes 

B 


VD 

—  >  25 
v 


(15) 


It  should  be  noted  that  this  order-of-magnitude  estimate  of  the  transition 
Reynolds  number  agrees  very  well  with  experimental  observations.  For 
example,  Viilu  [15]  found  the  value  of  11.2  tor  the  Reynolds  number  of 
breakdown  of  the  steady  laminar  jet.  Viilu' s  observations  were  later 
confirmed  by  Reynolds  [2]  who  reported  the  range  10  <  Re  <  30  for  the 
transition  Reynolds  number. 

The  buckling  number  criterion  of  transition  to  turbulence.  Equation 
(13),  explains  also  why  the  observed  transition  Re  is  a  number  consider¬ 
ably  greater  than  unity.  For  a  dimensionless  group  to  truly  delineate  the 
transition  from  one  mechanism  to  another,  it  must  have  a  value  of  order 
one  which,  after  all,  reflects  the  balance  between  the  competing  mechanisms. 
The  Reynolds  number  is  not  the  correct  dimensionless  group  to  describe 
transition  to  turbulence.  The  Reynolds  number  is  the  experimental 

(measurable)  reflection  of  the  N  transition  criterion. 

a 

Concluding  Remarks 

This  paper  unveiled  a  theoretical  basis  for  predicting  the  large 
scale  structure  and  fluctuating  behavior  of  inviscid  jets.  The  analogy 
between  jet  envelopes  and  elastic  columns  in  ax  al  compression  showed 
that  the  natural  tendency  of  inviscid  jets  is  to  buckle  over  a  precise 
wavelength  which  scales  only  with  the  jet  diameter.  The  repeated  buckling 
and  breakup  of  the  jet  column  is  responsible  for  the  sinusoidal,  river- like, 
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path  of  turbulent  jets  and  also  for  their  natural  whiplash  motion.  The 
predicted  natural  frequency  of  the  jet  agrees  very  well  with  measurements 
from  harmonic  excitation  experiments.  Finally,  we  learned  that  the 
transition  from  laminar  to  turbulent  jet  flow  occurs  when  the  buckling 
and  breakup  time  is  shorter  than  the  time  of  viscous  diffusion  across 
the  jet.  This  last  conclusion  is  supported  strongly  by  experimental 
observations. 

The  buckling  property  desc  ribed  in  this  paper  emerges  as  the  funda¬ 
mental  property  serving  as  origin  for  turbulent  motion  in  inviscid  flows. 
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Paper  no.  2 


The  buckling  property  of  viscid  and  inviscid  fluid  layers:  * 

Adrian  Be j  an 

Department  of  Mechanical  Engineering 
University  of  Colorado,  Boulder,  Colorado  80309 

Suleiman  and  Munson^  recently  presented  a  sequence  of  interesting 

experiments  concerning  the  buckling  of  thin  viscous  fluid  layers.  The 

authors  referred  in  their  paper  to  Taylor's  observations  on  the  buckling 

2 

of  slender  viscous  filaments  .  Taylor's  early  observations  and  the 
Suleiman  and  Munson  experiments  point  towards  an  important  analogy 
between  the  buckling  of  elastic  columns  and  the  buckling  of  viscous 
filaments  in  longitudinal  compression. 

The  object  of  this  letter  is  two-fold.  First,  it  brings  to  the 
readers '  attention  the  fact  that  a  buckling  theory  of  thin  viscous 
layers  already  exists.  Second,  it  points  out  that  the  buckling  of 
slender  columns  in  axial  compression  is  not  a  property  of  only  elastic 
solids  and  highly  viscous  fluids,  but  also  a  property  of  inviscid 
columns  (streams) . 

The  buckling  of  a  ;hin  viscous  layer  was  considered  theoretically 
by  Buckmaster,  Nachman  and  Ting3.  They  showed  that  viscous  layers  in 
longitudinal  compression  satisfy  two  basic  requirements: 

(i)  the  viscous  layer  is  in  axial  compression  and  the  compressive 
force  is  proportional  to  the  relative  velocity  between  the 
two  ends  of  the  layer; 

(ii)  if  curved,  the  viscous  layer  develops  in  its  cross-section  a 
bending  mcment  which  is  proportional  to  the  time-rate  of 
change  in  the  local  curvature. 


^published  In  Physics  of  Fluids.  Vol.  24,  pp.  1764,  1765,  1981. 
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The  momentum  equations  integrated  over  the  viscous  layer  led  Buckmaster, 

et  al. 3  to  a  global  equation  for  the  evolution  of  the  layer  centerline. 

This  equation  was  solved  by  assuming  various  initial  disturbances  as 

the  starting  shape  in  the  evolution  of  the  viscous  layer.  In  a  subse- 

4 

quent  paper,  Buckmaster  and  Nachman  have  extended  this  theory  to  the 
case  where  surface  tension  effects  play  an  important  role. 

Of  fundamental  interest  are  the  similarities  between  the  bucklinq 
of  a  viscous  layer  (the  viscida  problem)  and  the  buckling  of  a  slender 
elastic  column  (the  elastica  problem) .  A  close  examination  of  the 
Euler  theory  for  slender  elastic  columns  reveals  that,  for  infinitesimal 
buckling,  elastic  columns  must  satisfy  only  two  conditions5: 

(iii)  the  column  is  straight  and  in  axial  compression; 

(iv)  if  subjected  to  a  separate  bending  test,  the  column  develops 
in  its  cross-section  a  resistive  bending  moment  which  is 
proportional  to  the  local  curvature. 

The  infinitesimal  buckling  of  the  straight  column  follows  from  invoking 
static  equilibrium  for  the  system  sketched  in  Fig.  1.  For  a  straight 
or  nearly  straight  column  there  are  two  equilibrium  conditions  to  con¬ 
sider.  The  first  is  the  obvious  balance  of  compressive  forces  C  in 
the  longitudinal  direction.  The  second  condition  is  one  of  rotational 
equilibrium  which  states  that  the  excentricity  bending  moment  Cy  — 
however  small  --  must  be  balanced  at  all  times  by  the  cross-sectional 
bending  moment  M.  In  the  case  of  a  column  containing  elastic  material, 
a  special  bending  test  of  proscribed  curvature^  combined  with  knowledge 
of  the  elastic  properties  of  individual  fibers  in  the  slender  column 
leads  to  the  notion  of  a  cross-section  bending  moment  proportional  to 
the  local  curvature  (iv) .  Combining  this  notion  with  the  rotational 
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equilibrium  condition  yields  the  equilibrium  shape  of  the  straight 
column,  namely,  a  sinusoid  of  infinitely  small  amplitude.  The  infinitesi¬ 
mal  amplitude  is  indeterminate,  hence,  the  equilibrium  of  the  nearly 
straight  elastic  column  is  indifferent.  If  present,  the  slightest 
lateral  force  is  able  to  push  the  column  away  from  the  straight  equili¬ 
brium  shape  into  a  nearly  straight  (sinusoidal)  equilibrium  shape5. 

It  is  important  to  include  in  this  discussion  the  case  of  inviscid 
fluid  columns:  do  such  columns  conform  to  buckling  of  type  (i,  ii)  or 
to  buckling  of  type  (iii,  iv) ?  Thermodynamic  reasoning  alone  suggests 
that  inviscid  columns  should  buckle  according  to  model  (iii,  iv)  because, 

like  elastic  solids  and  unlike  highly  viscous  fluids,  inviscid  fluids 

7 

are  free  of  entropy  generation  . 

Referring  again  to  Fig.  1,  we  note  that  a  column  containing  inviscid 
fluid  is  none  other  them  the  imaginary  control  surface  drawn  around  an 
inviscid  jet  flowing  through  an  inviscid  fluid  at  rest.  Since  the 
column  is  stationary,  its  equilibrium  is  described  by  the  two  static 
conditions  (force  and  moment)  discussed  in  connection  with  the  Euler 
buckling  of  elastic  columns.  The  remaining  problem  is  to  determine  what 
special  forms  C  and  M  take  in  the  case  of  inviscid  stream  columns. 

For  the  compressive  force  on  the  control  volume,  C  ,  we  know  from  the 

Q 

thermodynamics  of  open  flow  systems  that 

C  =  p  A  V2  (1) 

where  p  ,  A  ,  V  are  the  density,  cross-sectional  area  and  velocity  of 
the  stream.  In  order  to  determine  the  cross-sectional  bending  moment 
M  we  conduct  a  special  bending  experiment  where  the  inviscid  stream  is 
held  in  a  duct  of  known  radius  of  cm^vature,  Rm  .  In  the  limit  of 


23 


vanishingly  small  curvature,  the  inviscid  flow  equations  show  that  the 

9 

resistive  bending  moment  in  the  cross-section  is 

M  *  p  V^I/R^  (2) 

where  I  is  the  area  moment  of  inertia  of  the  cross-section  . 

Equations  (1,2)  show  that  inviscid  fluid  columns  obey  conditions  (iii,iv) 
for  infinitesimal  buckling.  The  corresponding  equilibrium  shape  of  the 

9 

column  is  a  sinusoid  of  infinitely  small  amplitude  ,  whose  wavelength 
is 


X  =  2rr/l/A  (3) 

Inviscid  streams,  like  elastic  rods  and  viscous  layers,  possess  the 

natural  property  of  buckling.  The  wavelength  of  the  buckled  shape  scales 

only  with  the  transversal  dimension  of  stream:  for  example,  equation 

(3)  for  a  jet  of  round  cross-section  yields  X/D  «*  ir/2  .  Unlike  buckled 

5 

elastic  columns,  whose  equilibrium  is  indifferent  ,  buckled  inviscid 

columns  are  unstable.  The  post-buckling  evolution  of  the  inviscid  column 

is  well  understood,  forming  the  subject,  of  the  classical  theory  of  hydro- 
10 

dynamic  stability  .  However,  the  buckling  property  is  to  be  recognized 
as  responsible  for  the  wave-like  "disturbance"  assumed  routinely 
(empirically)  as  starting  point  in  any  hydrodynamic  stability  analysis. 

In  essence,  the  buckling  property  of  inviscid  streams  guarantees 
that  such  streams  cannot  flow  straight  through  another  fluid  or  through 
a  flexible  duct.  This  new  property  of  inviscid  fluids  stands  at  the 
very  root  of  the  phenomenon  of  turbulence.  For  example,  the  interaction 
between  the  elbows  of  the  buckled  stream  and  the  stagnant  medium,  leads 
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the  periodic  formation  of  large  eddies.  These  large-scale  turbulent 

structures  are  responsible  for  the  river-like  shape  of  turbulent  jets, 

wakes  and  plumes:  it  is  widely  observed  that  the  wavelength  of  this 

large-scale  meandering  path  scales  with  the  diameter  of  the  stream,  as 

predicted  in  equation  (3).  The  natural  buckling  of  inviscid  streams  is 

particularly  visible  in  the  case  of  streams  flowing  within  flexible 

boundaries,  like  rivers  and  capillary  rivulets  at  high  Reynolds  numbers. 

The  buckling  property  serves  as  theoretical  basis  for  predicting  other 

turbulence  parameters,  for  example,  the  natural  frequency  of  turbulent 

jets  (the  Strouhal  number)  and  the  critical  Reynolds  number  for  transi- 
9 

tion  to  turbulence  . 
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A.  Bejar.,  Buckling  and  Vibration  Theory  of  Jets,  Report  CUMER 
30-4,  Department  of  Mechanical  Engineering,  University  of  Colorado, 
Boulder,  1980. 

10  L.  Prandtl,  Essentials  of  Fluid  Dynamics  (Blackie  &  Son,  London, 
1969)  ,  p.  51. 


(control  volume)  In  axial  compression,  (b)  The  natural  buckling  wavelength 
of  circular  and  two-dimensional  inviscld  jets. 
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Paper  no.  3 


THEORY  OF  INSTANTANEOUS  SINUOUS  STRUCTURE 
IN  TURBULENT  BUOYANT  PLUMES  * 


by 


Adrian  Bejan 


Department  of  Mechanical  Engineering 
University  of  Colorado 
Boulder,  Colorado  80309 


This  paper  reports  a  theory  which  explains  the  flickering  motion 
of  turbulent  plumes  as  well  as  their  large-scale  sinuous  structure. 

The  theory  is  based  on  the  fact  that  the  invi scid  region  of  the  plume 
(the  "plume  column")  possesses  elastic  properties  analogous  to  those  of 
elastic  rods  subjected  to  longitudinal  oxnpression.  It  is  shown  that 
the  straight  plume  column  is  not  stable  and  buckles.  The  distance  between 
two  consecutive  elbows  is  proportional  to  the  local  plume  diameter,  in 
other  words,  the  shapes  of  all  buckled  plumes  are  geometrically  similar. 

A  buckled  plume  collapses  periodically  due  to  the  interaction  of  its 
lateral  elbows  with  the  stagnant  ambient.  This  interaction  is  responsi¬ 
ble  for  the  intermittent  formation  of  large-scale  buoyant  eddies  on  the 
periphery  of  the  turbulent  plume. 


* 

to  appear  In  Warme-und  Stof fubertragung,  Vol.  16,  1982. 
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1 .  Introduction 

Turbulent  plumes  represent  an  extremely  frequent  natural  phenomenon, 
in  the  atmosphere  as  well  as  in  the  hydrosphere.  Furthermore,  in  heat 
transfer  and  environmental  engineering,  the  turbulent  plume  constitutes 
one  of  the  most  effective  mixing  mechanisms  known  to  man.  The  state  of 
knowledge  on  turbulent  plumes  has  been  summarized  by  Turner  [1]  and,  as 
part  of  a  comprehensive  monograph  on  buoyancy-driven  flows,  by  Turner  [2]. 
Another  extensive  overview  of  the  field  was  published  most  recently  by 
Fischer,  List,  Koh,  Imberger  and  Brooks  [3].  A  critical  examination  of 
these  review  works  leads  to  the  conclusion  that  the  favored  approach  in 
turbulent  plume  research  is  one  where  the  "complications"  of  the  turbu¬ 
lent  flow  field  can  be  smoothed  out  by  the  classic  method  of  time¬ 
averaging  (Townsend,  (4]).  Consequently,  most  of  the  theoretical  work 
on  buoyant  jets  relies  on  the  axisymmetric  flow  model  which  lends  itself 
very  well  to  boundary  layer-type  analysis. 

Although  mathematically  attractive,  the  cime-averaged  turbulent 
t  i ume  concept  does  gross  injustice  to  the  physics  of  the  phenomenon.  A 
turbulent  exhaust  plume  does  not  rise  straight  up  into  the  air,  with  a 
fixed  shape  resembling  that  of  an  inverted  cone  or  funnel.  A  real 
plume  executes  a  periodic  lateral  movement;  the  plume  is  not  straight, 
rather,  it  has  recognizable  bends  separated  by  distances  of  the  ^ame 
order  as  the  plume  diameter.  The  smoke  rising  from  a  barbecue  or  camp 
fire  "snakes"  into  the  air.  This  large  scale  motion,  which  is  so 
effectively  blacked  out  by  time-averaging,  can  be  seen  clearly  in  Fig.  1. 

The  object  of  this  paper  is  to  present  a  theory  which  explains 
for  the  first  time  the  instantaneous  structure  and  the  "flickering" 
nature  of  turbulent  plumes. 


2.  The  Concept  of  Elastic  Plume  Column 


Consider  a  fluid  of  density  p  rising  vertically  through  a  heavier 
fluid  of  density  p  +  Ap  .  We  model  the  starting  section  of  this  buoyant 
stream  as  inviscid.  This  means  that  for  a  certain  height  at  the  base 
of  the  plume  we  are  neglecting  the  shear  interaction  between  the  plume 
fluid  and  the  stagnant  ambient.  This  is  permissible  for  a  height  of  the 
order  of  4-6  times  the  plume  diameter,  where,  as  shown  by  Crow  and 
Champagne  [5],  the  plume-ambient  shear  layer  is  chin  compared  with  the 
plume  diameter.  Having  neglected  the  interaction  at  the  plume-ambient 
interface,  it  is  helpful  to  imagine  the  plume  fluid  as  being  surrounded 
by  a  thin  flexible  sleeve  (control  surface)  which  is  stationary  with 
respect  to  the  stagnant  ambient. 

Imagine  now  a  cut  perpendicular  to  the  plume  axis.  If  the  column 
is  straight,  in  the  cross-section  generated  by  the  cut  we  distinguish  only, 
the  compressive  force  (impulse  [6]) 

C  =  ||  PV2  dA  ,  (2.1) 

A 

where  V  and  A  are  the  plume  velocity  and  cross-sectional  area.  If 
the  plume  column  is  slightly  curved,  we  distinguish  also  a  cross- 
sectional  bending  moment  M  .  This  banding  moment  is  due  primarily  to 
the  fact  that  the  plume  fluid  travels  relatively  faster  through  regions 
of  the  cut  which  are  located  on  the  inside  bank  of  the  bend.  It  is 
easy  to  show  that  the  cross-sectional  bending  moment  M  is  present  as 
soon  as  the  column  acquires  curvature.  Consider  for  this  purpose  a 
bending  test  in  which  the  uniform  stream  (p,  V,  A)  passes  through  a 
slightly  curved  ducts  the  radius  of  curvature  of  the  duct  centerline 
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,  is  infinitely  greater  than  the  plume  diameter,  D  .  The  Bernoulli 
equation  applied  along  a  streamline  requires  (Prandtl  [7]), 


+  P 


(2.2) 


where  the  left  hand  side  applies  to  the  straight  section  (P  =  0)  while 
the  right  hand  side  corresponds  to  the  curved  section.  Equilibrium  in 
the  radial  direction  requires  also  (Prandtl  (8] ) , 


dP 

dy 


(2.3) 


where  y  is  the  radial  coordinate  measured  in  the  plane  of  curvature 
away  from  the  plume  centerline,  toward  the  center  of  curvature.  Combining 
Equation  (2.2)  and  (2.3),  we  obtain 

v  =  V(1  +  -*-)  (2.4) 


The  cross-sectional  bending  moment  M  can  now  be  calculated  by  writing 


M  = 


(pv  +  P)  y  dA 


(2.5) 


which,  using  Equations  (2.4)  and  (2.2),  yields 


M  = 


PV 


z  r 


y  dA  . 


(2.6) 


Note  that  the  area  integral  appearing  in  Equation  (2.6)  is  the  area  moment 
of  inertia  (I  =  ttd4/64)  employed  routinely  in  the  study  o;  elastic  beam 
flexure  (Den  Hartog  [9]).  The  group  pv2  plays  the  role  of  "modulus  of 
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elasticity":  this  is  why  it  becomes  increasingly  difficult  to  manually 

bend  a  hose  as  the  turbulent  (inviscid)  flow  rate  through  it  increases. 

At  this  point  we  conclude  that  the  plume  column  is  in  axial  com¬ 
pression  and  that  the  cross-sectional  bending  moment  is  proportional  to 
the  curvature  of  the  centerline  (1/R^)  .  With  such  prepe’.ties,  the  straight 
plume  column  becomes  analogous  to  an  elastic  column  subjected  to  axial 
compression  (Den  Hartog  [10]).  This  analogy  guarantees  that  if  the  plume 
is  tall  enougn,  it  will  buckle  like  an  elastic  rod  compressed  between  the 
ends.  In  the  next  section  we  develop  the  "buckling  theory"  of  plumes 
viewed  as  elastic  columns. 

3.  The  Natural  Buckling  of  a  Straight  Plume  Rising  from  Rest 

Consider  the  stationary  control  volume  defined  by  two  cuts  normal 
to  the  plume  column,  at  z  and  z  +  dz  (Fig.  2).  The  balance  of 
vertical  forces  acting  on  this  volume  dictates 

^=ApAg  .  (3-D 

dz 

The  conservation  of  mass  through  the  control  volume  requires  also 

V (z) A (z)  =  Q  ,  constant  ,  (3.2) 

where  Q  is  the  volumetric  flowrate  through  the  plume  column.  Com¬ 
bining  these  conservation  statements  we  obtain 

V(z)  *  (2  ^  gz)**  ,  A (z)  =  Q/V(z)  ,  (3-3) 

where  we  made  the  additional  assumption  that  the  plume  fluid  originates 


from  re9t,  V(0)  =  0  . 
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As  in  the  buckling  theory  of  a  vertical  flagpole  (Timoshenko  and 
Gere  (11];  Den  Hartog  [12]),  the  rotational  equilibrium  of  the  plume 
element  dz  requires 

CdY  +  M  -  (M  +  dM)  »  0  (3.4) 


or,  using  Equations  (2.1)  and  (2.6), 


eAv2  S  +  ST  (pv2j  ==*»  =  0 


dz2 


(3.5) 


2  ** 

Note  that  in  the  limit  of  vanishingly  small  curvature,  1/R^,  =  -d  Y/dz“  . 
Taking  into  account  the  z-dependence  of  V  and  A,  Equations  (3.3), 
the  equati  n  of  flexure  (3.5)  becomes 


d3Y  4ir(2  g  Ap/p)**  h  dY 
dz3  Q  dz 


(3.6) 


The  first  important  conclusion  of  this  analysis  is  that  the  plume 
column  develops  a  "structure"  whose  characteristic  vertical  dimension 
scales  with 


z 

o 


2 

g  Ap/pJ 


1/5 


(3.7) 


The  vertical  length  scale  can  also  be  written  as 


z 

o 


zd“  1/5 
v  256' 


constant 


(3.7a) 


where  D(z)  is  the  plume  diameter  at  level  z  (note  that  D(z)  ~  z  )  . 
Introducing  the  dimensionless  coordinate  x  =  z/z  ,  the  equation  of 


flexure  becomes 


r-'*,  %■  w"1 
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MI 


d3Y  x  l]dy 
dz3  X  dx 


0 


(3.8) 


The  general  solution  to  this  equation  is  expressible  in  terms  of  Bessel 
functions  (Watson  [13]) 


£  -  *,  *1/2  Ja/5(4/5  x5/4)  ♦  K2  x1/2  J_2/5(4/5  «V4, 


(3.9) 


where  and  K2  are  arbitrary  constants.  Noting  that 

lim  =  Kxx  +  k2  *  (3.10) 

x-*-0 


we  set  K2  =  0  ,  which  accounts  for  the  fact  that  in  the  very  beginning 
the  buoyant  fluid  rises  vertically.  The  shape  of  the  plume  column,  Y(x)  , 
is  obtained  by  integrating  (3.9), 

? 

1  y5 

Y  (x)  =  K  m  '  J„._(m)  dm  (3.11) 

2/5 

0 


where 

£  =  4/5  x5/4  .  (3.12) 

The  integral  appearing  in  (3.11)  was  calculated  with  sixth-digit  accuracy, 
by  first  expanding  J  (m)  in  a  power  series.  The  principal  values  of 
this  integral  are  reported  here  in  Table  1,  because  they  are  not  available 
in  the  literature. 

The  magnified  shape  of  the  buckled  plume  column  is  shown  in  Fig.  3 
as  the  function  Y(x)/K.  It  should  be  remembered  that  Equation  (3.11)  is 
the  result  of  a  small  amplitude  analysis  (|dY/dz|  «  1)  ,  hence,  the 

amplitude  K  is  negligible  when  compared  with  the  vertical  length  scale 
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z  .  The  static  equilibrium  of  the  straight  plume  column  is  indifferent 
o 

[ 10 1 ;  however,  given  the  slightest  lateral  disturbance,  the  plume  buckles 
in  accordance  with  t.ie  vertical  periodicity  shown  in  Fig.  3.  The  buckled 
shape  is  a  succession  of  elbows  whose  spatial  frequency  increases  with 
height.  The  first  four  elbows  are  located  at 

x.  =  2.878,  x„  *  5.101,  x,  =  7.101,  x„  «  8.967  . 

12  3  4 

According  to  the  classical  argumr.nts  of  Hydrodynamic  Stability  Theory 
[14],  the  sinuous  shape  will  be  amplified  due  to  the  dynamic  interaction 
of  the  elbows  with  the  stagnant  ambient.  It  is  this  interaction  which 
leads  to  the  formation  of  large  turbulent  eddies  visible  in  most  atmospheric 
plumes  (Fig.  1) . 

4.  Formation  and  Evolution  of  Large  Eddies 

The  mechanism  of  elbow  turbulence  f <  ;ation  is  shown  schematically 
in  Fig.  4.  Speaking  only  qualitatively,  in  the  first  phase  of  the  pro¬ 
cess  the  segment  of  plume  column  located  in  the  vicinity  of  a  natural 
ellxiw  is  pushed  to  the  side  by  the  horizontal  resultant  of  the  two  axial 
corr.oressive  forces.  The  lateral  movement  of  the  elbow  is  eventually 
blocked  by  the  pressure  build-up  associated  with  either  the  stagnation  of 
some  plume  fluid  into  the  still  ambient,  or  with  accelerating  some  of  the 
ambient  fluid.  This  pressure  build-up  is  consistent  with  the  fact  that 
in  a  column  of  finite  curvature  the  pressure  increases  in  the  radial 
direction  (see  Equation  (2.3)].  AC  the  end  of  the  first  phase  of  the 
process,  the  elbow  has  degenerated  into  a  buoyant  eddy  of  mixed  fluid? 
this  eddy  rotates  as  shown  in  the  figure,  and  rises  slower  than  the 


.  v:  • 
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unmixed  plume  fluid  found  at  the  same  altitude. 

The  second  phase  of  the  mechanism  is  triggered  by  the  effect  of 
high  pressure  nodules  already  present  on  both  sides  of  the  plume  column: 
the  plume  buckles  in  a  new  mode,  feeding  its  stream  through  the  struc¬ 
tures  generated  by  the  old  set  of  elbows.  At  the  end  of  the  second 
phase,  the  new  column  gives  birth  to  a  new  set  of  large  eddies  which  con¬ 
tinue  to  rise.  In  general,  the  eddies  produced  by  the  upper  elbows  inter¬ 
fere  with  older  eddies  rising  from  lower  levels.  This  interaction  is 
partly  responsible  for  lateral  growth  of  the  turbulent  plume,  as  x 
increases . 

The  most  regular  feature  in  the  evolution  of  the  turbulent  plume 
column  is  the  root  section,  0  <  x  <  ,  which  precedes  the  formation  of 

the  first  elbow.  The  root  of  the  plume  swings  back  and  forth,  not  neces¬ 
sarily  in  che  same  plane,  with  a  characteristic  frequency,  f^  =  l/(2At) . 
The  half  period  At  of  this  flickering  motion  can  be  calculated  by  inte¬ 
grating  the  first  of  equations  (3.3):  At  is  the  time  interval  needed 
by  plume  fluid  co  rise  from  z=0  to  the  first  elbow,  z,  =  x, z  . 
Combining  this  result  with  Equations  (3.7)  and  (3.3)  ,  we  discover  that 
the  Strouhal  number  based  on  the  plume  velocity  and  diameter  at  the  first 
elbow,  V  and  ,  is  a  universal  constant  for  all  inviscid  plumes 

f  1D1 

St  =  =  0.534  .  (4.1) 

1  V1 

Another  interesting  result  of  the  buckling  theory  is  that  the  large 
elbow  structure  (the  distance  between  elbows)  scales  only  with  the  plume 
diameter,  which  means  that  the  buckled  shapes  of  all  inviscid  plumes  are 
geometrically  similar.  This  prediction  agrees  with  the  observed  structure 
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of  atmospheric  plumes.  According  to  Equation  (3.7a) ,  the  plume  diameters 
corresponding  to  the  first  four  elbows  are, 

—  -  =  3.071,  2.662,  2.45,  2.332  . 

z  % 

O  X. 

1 

5  shows  a  scale  drawing  of  a  round  plume  and  the  relative  location 
first  four  elbows.  This  universal  geometry  agrees  very  well  with 
observations  of  turbulent  plumes  rising  above  barbecues  and  camp 
It  also  agrees  with  the  plume  photographed  in  Fig.  1. 

5.  Laminar  Plumes  vs.  Turbulent  Plumes 

The  instantaneous  sinuous  structure  and  the  eddy  formation  mechanism 
described  so  far,  apply  only  to  large  Reynolds  number  (turbulent)  plumes 
which  can  be  modeled  as  inviscid.  The  present  theory  does  not  apply  to 
laminar  plumes.  It  is  important  to  note,  however,  that  the  universality 
of  the  sinuous  structure  (Fig.  5)  provides  a  theoretical  explanation 
for  the  phenomenon  of  transition  from  laminar  to  turbulent  plume  flow. 

According  to  Equation  (4.1),  the  base  of  the  plume  column  fluctuates 
within  a  characteristic  time  interval  which  is  proportional  to  the  dia¬ 
meter  and  inversely  proportional  to  the  velocity.  This  local  characteristic 
time, 


Figure 
of  the 
visual 
fires. 


t 


fluctuation 


D  /V 

r  1 

St, 


(5.1) 


represents  the  "natural  heartbeat"  of  the  plume  as  an  inviscid  stream. 
If  the  plume  flow  is  highly  viscous  (laminar) ,  the  plume  column  has  a 


-37- 


different  characteristic  time:  the  viscous  local  time  is  associated 
with  the  travel  of  information  by  viscous  diffusion  from  the  plume-ambient 
interface  to  the  plume  centerline  (total  distance  D/?. )  .  The  viscous 
communication  time  is  given  by  the  solution  to  Stokes'  first  problem  [15], 


V2 


2  (vt  .  ) 

viscous 


1  or  t  . 

viscous 


16v 


(5.2) 


The  plume  flow  remains  laminar  if  the  ambient  can  communicate  viscously 
with  the  stream  faster  than  the  stream  can  fluctuate.  If  the  column  can 
fluctuate  faster  than  the  viscous  diffusion  time,  then,  of  course,  the 
flow  opts  for  the  inviscid  (turbulent)  regime.  At  transition,  the  viscid 
and  inviscid  time  scales  are  of  the  same  order  of  magnitude;  combining 
Equations  (5.1)  and  (5.2)  we  conclude  that  the  laminar-turbulent  transi¬ 
tion  is  marked  by 


V  D 
11 


30 


(5.3) 


This  theoretical  estimate  of  the  transition  Reynolds  number  agrees 
very  well  with  experimental  observations  and  with  estimates  based  on 
hydrodynamic  stability  analysis.  Experimentally,  it  is  well  known  that 
the  transition  in  free  jet  and  wake  flow  occurs  in  the  vicinity  of 
Re  -  30  (see,  for  example,  Schlichting  [16]).  The  hydrodynamic  stability 
of  buoyant  plumes  and  wall  layers  has  been  studied  extensively,  as 
summarized  by  Gebhart  [  17 ] .  For  example,  the  lowest  Reynolds  number 
where  instability  has  been  detected  in  vertical  laminar  boundary  layer 
flow  over  a  constant-flux  wall  is  of  order  67  (see  Fig.  8-21  in 
Ref.  [ 17 ] ) .  The  difference  between  this  estimate  and  Equation  (5.3)  is 
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explained  by  the  fact  that  in  Ref.  [ 1 7 ]  the  local  Reynolds  number  (G*) 
is  based  on  the  thickness  of  the  laminar  boundary  layer.  If,  instead, 
the  Reynolds  number  is  expressed  based  on  the  displacement  or  momentum 
thickness  (more  consistent  with  the  slug  flow  model  employed  in  the  pre¬ 
sent  theory)  then  the  transition  Re  is  a  number  of  the  same  order  as  in 
Equation  (5.3). 


6.  Concluding  Remarks 

The  demonstrated  ability  to  correctly  predict  the  transition  to 
turbulence  is  additional  evidence  that  the  sinuous  structure  (Fig.  5) 
and  the  characteristic  time  scale  (Eq.  4.1)  are  real  properties  of  all 
turbulent  buoyant  plumes.  It  is  important  to  recognize  also  that  the 
large-scale  buckling  phenomenon  described  in  this  paper  is  not  observed 
only  in  buoyant  pxumes  as  in  Fig.  1.  The  same  phenomenon  has  been 
observed  and  described  as  "orderly  structure",  "corkscrew  shape",  and 
"whiplash  motion"  in  the  starting  section  of  low  speed  turbulent  jets 
(Crow  and  Champagne  [5];  Reynolds  [ 18 ] ) .  Others  have  recognized  this 
snake-like  shape  as  the  cause  of  the  intermittency  phenomenon  (Yih  [ 19 ] , 
Fig.  8b,  and  pp.  545,  546).  In  the  field  of  boiling  heat  transfer  we  are 
familiar  with  the  occurrence  of  tall  S-shaped  vapor  bubbles  (buckled 
vapor  plumes)  on  intensely  heated  horizontal  surfaces:  this  phenomenon  is 
recognized  as  the  "continuous  vapor  column  regime"  (see  Figs.  2(d)  and 
2 ( j )  in  Moissis  and  Berenson  [20]). 

The  reader  can  easily  reproduce  the  natural  buckled  shape  by  experi¬ 
menting  with  the  continuous  water  column  falling  from  the  kitchen  faucet. 
This  water  column  is  a  "sinking"  plume  defined  not  by  a  flexible  plume- 
ambient  interface  as  in  Fig.  1,  but  by  the  flexible  hose  provided  by 
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capillary  forces.  Placing  his  finger  about  1-2  cm  under  the  faucet,  the 
reader  can  buckle  this  sinking  plume  into  a  shape  which  resembles  very 
closely  the  shape  shown  in  Fig.  5,  rotated  by  180°.  tn  the  kitd.en 
faucet  experiment  the  buckled  column  does  not  break  up,  since  the  lateral 
growth  of  its  elbows  is  suppressed  by  the  effect  of  surface  tension. 

The  same  stabilizing  effect  is  present  in  the  S-shaped  tall  bubbles 
photographed  during  intense  boiling  by  Moissis  and  Berenson  [20 ] . 

From  the  point  of  view  of  theoretical  research  in  fluid  mechanics, 
the  Buckling  Theory  reported  in  this  paper  represents  a  substitute  for 
Hydrodynamic  Stability  Theory  [ 1 7 ] ,  for  determining  the  instability  of 
large  Reynolds  number  plume  flow.  This  conclusion  is  not  surprising, 
since  the  Navier-Stokes  equations  used  in  Hydrodynamic  Stability  Theory, 
and  Equations  (2.3),  (3.1)  and  (3.4)  used  in  the  Buckling  Tneory,  are 
statements  of  the  same  (classical)  principles  of  mechanical  equilibrium. 
The  transition  to  turbulent  plume  flow,  discussed  in  Section  5,  indicates 
that  the  natural  buckled  shape  identified  by  the  Buckling  Theory  corres¬ 
ponds  to  the  least  stable  asymmetric  disturbance  determined  via  Hydro- 
dynamic  Stability  Theory. 
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Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


List  of  Captions 

Large-scale  periodic  structure  in  the  plume  above  a  natural 
gas  well  on  fire  —  Douglas  Pass,  Colorado  (AP  Laserphoto, 
reprinted  from  the  Longmont  Daily  Times-Call,  Oct.  18-19,  1980) . 

Static  equilibrium  of  a  segment  of  vertical  plume  column. 

The  buckled  shape  of  the  plume  column  centerline. 

Mechanism  of  formation  of  large-scale  eddies  in  a  turbulent 
plume  rising  from  rest. 


Figure  5. 


Universal  geometry  of  buckled  plume  columns. 
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